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Abstract 

We study the behavior of the two-dimensional two-component plasma 
in the presence of some adsorbing impurities. Using a solvable model, we 
find analytic expressions for the thermodynamic properties of the plasma 
such as the n-body densities, the grand potential, and the pressure. We 
specialize in the case where there are one or two adsorbing point impurities 
in the plasma, and in the case where there are one or two parallel adsorbing 
lines. In the former case we study the effective interaction between the 
impurities, due to the charge redistribution around them. The latter case 
is a model for electrodes with adsorbing sticky sites on their surface. 
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1 Introduction 

A two-component plasma (TCP) is a system composed of two different types 
of particles with charges ±e interacting through a Coulomb potential. In this 
work, we study a two-dimensional model which is exactly solvable. The solution 
of Poisson equation shows that in two dimensions, the Coulomb interaction 
between two particles, with charges q±, q 2 , at a distance r from each other, is 

v(r) = - qi q 2 ln(r/L), (1.1) 

where L is an irrelevant length scale fixing the potential reference. 

Let us define the dimensionless coupling constant T — (3e 2 , where f3 — 
(fcgT) -1 , T the temperature and fcs the Boltzmann constant. The coupling 
constant T is the ratio between the electrostatic and thermal energy of the 
particles of the gas. If T < 2, the thermal energy is high enough to prevent 
the collapse of particles of different sign. If T > 2, the system is unstable 
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against the collapse and the point-particles must be replaced by hard disks of 
diameter a. Using the analogy between the two-component plasma and the sine- 
Gordon quantum field theory, exact results for the thermodynamic properties 
of the Coulomb gas, in the bulk [2] and near a planar interface [3, 4], have been 
obtained when T < 2. For r = 2, the classical two-dimensional two-component 
plasma is equivalent to a quantum free Fermi gas [1]. This fact allows us to find 
analytic expressions for the thermodynamic properties and correlation functions 
of the plasma, for a wider variety of geometries [1, 5, 6], and even in the presence 
of external fields [7]. In this work, we restrict our analysis to the case T = 2. 

Assuming T = 2, we pretend to find analytical solutions for the thermody- 
namics of the two-component plasma in the presence of some impurities. Results 
have been obtained for a similar problem, the one-component plasma with ad- 
sorbing impurities [8, 9]. For the two-component plasma, some similar problems 
to the present one, with adsorbing boundaries, were considered in Refs. [10, 11]. 

An impurity in our model can be understood as a particle, different from the 
charged particles of the gas. It interacts with the particles of the gas through 
a non-electrical potential U±(r). As usual, for inhomogeneous systems under 
the action of an external field derived from a potential U±(r), we can define a 
position depend fugacity due to the external potential created by the impurity 
m±(r) = me~ /3 ' 7± ( r ), where m is the bulk fugacity. Here we are interested in 
adsorbing impurities, which can be modeled with a position depend fugacity 
m±(r) = m[l + a±5(r — R)l , where a± is the magnitude of the interaction of 
the impurity with the positive and the negative particles of the plasma (it will 
be called adhesivity), while R is the position of the impurity in the plasma. In 
order to prevent the collapse of the opposite charged particles, we will consider 
that the impurity only attracts one type of particles of the system. For instance, 
if the impurity only attracts the positive particles of the system, a + = a and 
a_ = 0. 

This simplified model could be applied to some systems. An example is a 
salt solution with positive and negative ions and porous sites, which can attract 
some of the ions, with different sizes for each species. The size of the porous sites 
can be such that they are small enough in order to avoid that the larger ions 
get adsorbed into the sites, while the smaller ones are allowed to be adsorbed. 
The same salt solution may have some "dusty-points" ; the contact interaction 
may be generated by a chemical interaction among those points and the charged 
particles. An example of a chemical interaction can be presented in some atoms 
whose valence layer is not completely filled and they can receive electrons from 
an atom with a small electronegativity. 

The present document is organized as follows. In section 2, we briefly review 
the general mathematical treatment which allow us to find the thermodynamic 
properties and correlation functions of the two-component plasma in the pres- 
ence of any external potential when T = 2 [1]. In section 3, we consider the case 
where there are a few point impurities in the plasma. We concentrate specif- 
ically in the case of one or two impurities. We compute the density profiles 
around the impurities and the two-body correlation functions using the method 
of Ref. [1]. Also, we extend the treatment presented in Ref. [8] for the one- 
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component plasma, to find general expressions for the grand potential, and the 
one- and two-body density functions for the two-component plasma with impu- 
rities, in terms of the same quantities for the unperturbed system. In section 4, 
we study the two-component plasma with one or two lines of adsorbing impu- 
rities. We compute the one- and two-density functions and the grand potential 
using the method of Ref. [1]. A summary of the main results is presented in 
section 5. 



2 Method of solution 

To solve our specific model, we use the method introduced by Cornu and Jan- 
covici [1]. This mathematical treatment is valid for a globally neutral two- 
component Coulomb gas with a coupling constant T = 2. Since the model is a 
two dimensional one, it is convenient to express the position r = (x, y) of the 
particles in complex coordinates as z = x + iy. As explained earlier, at T = 2, 
it is mandatory to introduce a cutoff distance a between particles to avoid the 
collapse between particles of opposite sign. This cutoff can be understood as a 
small diameter of the particles. It is shown in Ref. [1] that, in the limit a — > 0, 
the grand partition function S in this scheme is given by 



InS = Tr 



In ( ^ + m + (r) 1 + <7z + ra_ (r) — 



(2.1) 



where = <J x d x + o y d y is the two dimensional Dirac operator. The matrices 
(J x ,a y ,a z are the Pauli matrices. The position dependent fugacity, defined as 
m±( r ) = me"'" ,± W, takes into account an external potential U±(r) which 
can act differently on the positive or negative particles. The constant rescaled 
fugacity m is related to the chemical potential [i by m = 2 r Ke^ ti Lj a 1 . It has 
units of inverse distance, and it is shown in Ref. [1], that mT 1 is the screening 
length. 

The calculation of the fc-body densities is reduced to find the set of Green 
functions G SlS2 (ri, r 2 ), with s±, s 2 = ±1, of the operator intervening in Eq. (2.1). 
They are the elements of a 2 x 2 matrix which satisfy the differential equations 

( m+(n) 2d z \f G ++ (r 1 ,r 2 ) G+_(ri,r 2 ) \ _ x( , . , 

The density n± (r) and two-body truncated density correlation function n± lT (r i , r 2 ) 
of the system are given by [1] 

n s (r) = m s (r)G ss (r, r) (2.3a) 
n^ 2 (ri,r 2 ) = -m Sl (r)m S2 (r 2 )G SlS2 (ri,r 2 )G S2Sl (r 2 ,ri), (2.3b) 

where T means truncated. More generally, the truncated fc-body densities are 
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given by [1] 

ni Ji )T s k ( r U ■ ■ ■ ,r k ) = {-l) k m Sl {v 1 )...m Sk {v k ) 



E 



G, 



(2.3c) 

( r U ! r i2 ) ' ' ' Gs ifc s il ( r i fe i r ii ) 



cycles (ii—ii,) 

The Green functions satisfy the useful relations 



G ss (ri,r 2 ) = G ss (r 2 ,ri) 
G s _ a (ri,r 2 ) = 



-G_ ss (r 2 ,ri). 
The grand potential (3Q = — In S can be written as 
/3fi = -^ln(l + A fc ) 



(2.4a) 
(2.4b) 



(2.5) 



where are the eigenvalues of the system 



m+(r) 
m_(r) 



ip(r) = 0, with 



^(r) 



_ f ff(r) 



(2.6) 



a two-component spinor. Equation (2.6) shows the mentioned equivalence be- 
tween the two component plasma at T — 2 and a free fermion gas. The pressure 
p is given as usual by 



<91nS 1 



(2.7) 



in the thermodynamic limit, where il is the grand potential. The "volume" of 
the gas is replaced by its area A because we are dealing with a two-dimensional 
model. 

The above results are valid in the continuum limit a — > 0. However, in 
this limit the partition function and other thcrmodynamical quantities such as 
the pressure and the density are divergent. We shall compute the (divergent) 
dominant term as a function of the cutoff a. On the other hand, the truncated 
correlation functions arc finite in the limit a — > 0. 



3 The plasma with some point adsorbing impu- 
rities 

In this section, we study the two-component plasma in the presence of one or two 
impurities. First, we find the Green functions associated to this system at T = 2. 
This allows us to find the fc-body densities of the plasma. We will distinguish 
two kind of impurities. They will be called "positive" impurities when they 
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attract the positive particles of the plasma, and "negative" impurities in the 
opposite case. This is a small abuse of language, since the impurities do not 
carry any electrical charge by themselves alone. 

Then we shall extend the general theory presented in Ref. [8] , to find general 
expressions for the thermodynamic properties of the system with impurities in 
term of those same quantities for the unperturbed system. These last results 
are valid for any value of the coupling constant T. 



3.1 Density and correlations at T = 2 
3.1.1 A positive impurity 

We consider a single impurity in the plasma that only adsorbs the positive 
charged particles. In this case the fugacity becomes m + (r) = m,[l+ai8(r— Ri)] 
and 77i_(r) = m. The constants a\ and Ri represent the adhcsivity of the 
impurity and its position. 

First we compute the Green functions that are needed to find the density 
and correlations. To solve the system (2.2) for the Green functions, we assume 
solutions of the form G SlS2 = G° 1S2 + Gl 1S2 , where G° 1S2 are the bulk solutions 
(solutions for an unperturbed plasma with m± (r) = m) [1] 

77? 

G° ± (n,r 2 ) = _Jf ( TO | ri -r 2 |), (3.1a) 

n a , \ m fa - x 2 ) + i(yi - y 2 ) „ , . 
G_ + (n,r 2 ) = ^— - { JTiMn-ral) 

= Je^^Hn-rjj), (3.1b) 

where #i 2 is the polar angle of the vector ri — r 2 and K (x), K\{x) are the 
modified Bessel functions of second kind of order and 1. The bulk solutions 
satisfy the appropriate boundary condition at ri = r 2 imposed by the 5{v\ — r 2 ) 
in Eq. (2.2). The functions G l SlS2 take into account the contribution of the 
impurity. They satisfy the differential equations 

mG^ + (r 1 ,r 2 ) + 2a zi Gi + (r 1 ,r 2 ) 

+™ ai <5(r 1 -R 1 )[G° + (r 1 ,r 2 )+G^ + (r 1 ,r 2 )] - (3.2a) 

2d 2l GV+(ri,r 2 ) + mGl + (n,r 2 ) = (3.2b) 

mG^_(r 1 ,r 2 ) + 29 Zl G__(r 1 ,r 2 ) 

+mai<y(ri-Ri)[G^_(ri,r 2 ) + G^_(ri,r 2 )] = (3.2c) 

2d- Zl G 1 + _(v 1 ,v 2 ) + mG 1 __(v 1 ,r 2 ) - 0. (3.2d) 

We will first solve the equations for G\ + and G 1 _ + . We divide the space into 
two regions: r\ < R\ and n > R\. These regions will be denoted by the 
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superscripts (1) and (2) respectively. The general solutions for G± + arc 



Gl + (n,r 2 ) 
Gi+(ri,r 2 ) 



Z l& e U9l 4PMmri), n > R 1 . 



n < Ri 



(3.3a) 
(3.3b) 



with Ii(x) and Ki(x) the modified Bessel functions of first and second kind of 
order I. From Eq. (3.2a), we notice that G 1 ++ is continuous, but that G 1 _ + must 
be discontinuous at ri = Ri, due to the <5(ri — Ri) term. Multiplying Eq. (3.2a) 
by e~ lWl and integrating on ri in a small annulus domain, centered at the origin 
and containing Ri, gives 



Gi+(ri,r 2 ) 



ri=R+ 



Gi+(ri,r 2 ) 



ri=R- 



£^[G° ++ (Ri,r 2 ) + G] h+ (R 1) r 2 )] = 

ri=RJ 



ri=R 1 



(3.4a) 



, (3.4b) 



where gL 1 + +1 refers to the I + 1 term of the sum in Eq. (3.3b). Imposing the 
boundary conditions (3.4), we find 



and 



4 (2) _ (i) Jj(mfii) 



TO 3 aiif;(mi?i)e- l(8 «i J Fi:o(™|Ri - r 2 |) 



4tt 2 [1 



■ K^mRxjl^mRx 



(3.5) 



(3.6) 



where 8^ is the polar angle of the vector Ri. 

At this point, let us comment a few details about the bulk density, for 
an unperturbed system, which will be useful to interpret the denominator in 
Eq. (3.6). The bulk density is given by n±(r) = mG° ±± (r, r). Nevertheless, 
this expression diverges, since, for small argument, Kq(x) ~ ln(2/x) — C, with 
C ~ 0.5772 the Euler constant [1]. This divergence can be avoided by replacing 
the point particles by hard spheres of size a, which represents the minimum 
distance between two charged particles. Then, we compute the bulk density as 
n± = raGj_ ± (ri, r 2 ) with |ri — r 2 | = a — > 0. Thus, 



n 



2tt 



K (m\r 1 - r 2 |) 



2tt 



Ko(ma) 



m 

2lT 



In- 



C 



ma 



(3.7) 



where we used the small- argument expansion of the Bessel function _Ko- Now, 
let us recall the expansion K (m\ri — r 2 |) = X)«ez e™^ 1 " 92 '/„(m7' < )]f n (mr > ) 
with r< = min(ri,r 2 ) and r> = max(ri,r 2 ) [12]. With this expansion, we can 
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obtain another formal (divergent) expression for the bulk density 



n° ± (Rx ) = n = — K n (mRi )I n {mR\ ) . 



(3.6 



Since the unperturbed system is homogeneous the bulk density is constant, 
n±(Ri) = no does not depend on Ri. We notice that it is this expression that 
precisely appears in the denominator of Eq. (3.6). Then, 



47r 2 (l + ain ) 



(3.9) 



Replacing in Eqs. (3.3a) and (3.3b), 



G\ + (r u r 2 ) = - 



m 3 ai 



47T 2 [l + OtlTlQ 



■if (m|Ri - r 2 |)i<ro(m|ri - Ri|) (3.10a) 



m 3 rtT p i(>1R i 

Gi + (n,r 2 ) = -—^ T K a (m\r 2 - R^^Hn - R^) (3.10b) 

A-k z [I + ain \ 

with 9ir ± the polar angle of the vector ri — Ri. 

In a similar way, the other two Green functions are 

Gi_(n,r 2 ) = \ 11 : — ifi(m|ri - Ri|)ifi(m|r 2 - Ri|) 



4tt 2 [l + awo] 



(3.11a) 



GV-(ri,r 2 ) = 



47T 2 [l + ain ] 



Jf (m|ri-Ri|)iri(m|Ri-r2|) (3.11b) 



with 92R ± the polar angle of the vector r 2 — Ri . The one-body densities can be 
calculated by using Eq. (2.3a), the result is 



n+(r) = [l + ai6(r - Ri) 



no 



_ /m 2 \ 2 Qi [^ (m|Ri-rl)] : 
V 2tt / [1 + am ] 



n_(r) 



"o + 



/m 2 \ 2 ai[^i(mjRi-r|)]' 



(3.12a) 
(3.12b) 



V 27r / [1 + ain J 

The number of positive and negative particles of the plasma, N± = J n± (r) d 2 r , 



is 



N+ = 7V_ = n nR 2 



m 2 a\ 



+ 



n ai 



47r[l + ain ] 1 + a i n o 



(3.13) 



where we supposed that the plasma is confined in a large disk of radius R. 

The difference N + — N- = 0. This shows that the global neutrality of the 
system is not changed, only the charge distribution. The amount of positive 
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Figure 1: Average number of adsorbed particles as a function of the adhcsivity. 
The cutoff is ma — 0.01. The dashed line is the asymptotic value = 1 for 
ai — > oo. 



charge bound to the impurity, TV" 1 , is given by the term which multiply the 
delta distribution in Eq. (3.12a), evaluated at Ri. The result is 

iV? 1 = _^2_, (3.14) 
+ 1 + amo 

which can be interpreted as a mean occupation number or probability that 
the adsorbing point is occupied by a positive charge. Figure 1 shows ./V" 1 as 
a function of ct\. As the adhesivity ct\ is increased, the average number of 
adsorbed positive particles A^" 1 increases, as expected. In the limit a\ — > oo, 
the maximum value of adsorbed particles is obtained N? 2 = 1. Notice that 
in average there cannot be more than one adsorbed particle per impurity site. 
Once a particle is adsorbed, any other particle of the same sign feels a strong 
electrostatic repulsion that prevent it from approaching the adsorbtion site. 

This localized positive adsorbed charge is the responsible of the change in 
the charge distribution. Since the bound charge is positive, it repels the positive 
particles and attracts the negative ones. As a consequence, the density of the 
negative particles around the impurity increases, and the density of the posi- 
tive particles around the impurity decreases. As it is expected, the larger the 
adhesivity ol\ is, the larger this effect is. This behavior can be seen in figures 2, 
3 and 4 which show the density profiles of the positive and negative particles, 
and the charge density profile. 

If the impurity attracts the negative particles, instead of the positive ones, 
we obtain the same results by interchanging (+ < — ► — ). 
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0.25 0.5 0.75 1 1.25 1^5 1.75 2 mr 

Figure 2: Density of positive particles as a function of the distance for m 2 a\ = 
1,0.5 and 0.2 from the darkest to the lightest. The dashed line represents the 
bulk density. The vertical line at mr — 1 represents the Dirac distribution 
related to the impurity position which is located at mRi = x. The angular 
difference is 9 — 9^ = 0. The cutoff is ma = 0.01. 
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0.25 



mr 



Figure 3: Density of negative particles as a function of the distance for m 2 a\ = 
1,0.5 and 0.2 from the darkest to the lightest. The dashed line represents the 
bulk density. The angular difference is 6 — 9^ = 0. The cutoff is ma = 0.01. 



m [n+ -n_ ] 
2r 



1 - 




Figure 4: Charge density profile as a function of the distance for m 2 ai = 1,0.5 
and 0.2 from the darkest to the lightest. The vertical line at mr — 1 repre- 
sents the Dirac distribution related to the impurity position which is located at 
mRi = x . The angular difference is 9 — 9^ = 0. The cutoff is ma = 0.01. 
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0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 5: n\ _ (ri, r2) as a function of r\ for m 2 ai = 10 and 0.5 from the darkest 
to lightest. We assumed mr 2 = 0.5 and mi?i = 1. The angular differences are 
&i — #2 = 0i — 0Ri = 0: we are looking through a line which passes over r 2 and 
Ri. The delta function at ri = Ri is not showed. 




-2.5; 

— — — — i — . — . — — . — . — . — — . — . — . — — . — — — TUT — _ — , — , — , — — , — mr 
0.2 0.4 0.6 0.8 1 1.2 1.4 

Figure 6: n^l T (ri,r 2 ) as a function of the distance. We assume the same 
parameter values as in figure 5. 
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m n 



(2) T 




0.2 0.4 0.6 0.8 



mr n 



Figure 7: n\ L (ri,r2) as a function of the distance. We assume the same 
parameter values as in figure 5. The delta function is not plotted. 



As seen in figures 2, 3 and 4, a charge polarization cloud is formed around 
the impurity. It is interesting to evaluate the electric potential created by this 
charge distribution. To compute the electric potential, we use the Poisson equa- 
tion V 2 0(r) = — 27re[n + (r) — n_(r)], where n±(r) are given in Eqs. (3.12). In 
the thermodynamic limit the plasma is large enough and we can assume an an- 
gular symmetry around the impurity (locating it at the origin). Thus, Poisson 
equation reduces to 

J|:[ r i>( r )] =-27re[n+(r)-n_(r)]. (3.15) 
Integrating, we find the electric potential 

*(r) = r e ™y T [Ko(mr)] 2 = 6 -^±- [K (mr)] 2 . (3.16) 
2ir[l + ain \ znon 

From Eq. (3.16) we can see that close to the impurity, the potential behaves as 
[ln(mr)] , which is stronger than the bare Coulomb potential — ln(r/L). This is 
due to the strong electrostatic coupling considered here [17], for smaller coupling, 
one would expect that the electric potential close to the impurity behaves as the 
bare Coulomb potential. Far away from it, the interaction decays exponentially 
due to the screening effect. 

Replacing the results for the Green functions, Eqs. (3.10) and (3.11), into 
Eq. (2.3b) allows us to obtain the two-body density correlation functions. Fig- 
ures 5, 6, and 7 show plots of the two-body densities, for different values of the 
adhesivity. 
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3.1.2 Two positive impurities 



Now, we consider two impurities located at Hi and R 2 with adhesivities ol\ and 
a 2 respectively, both attracting positive particles. The fugacities are now given 
by m + (r) = m[l + ai<5(r — Ri) + a 2 S(r — R 2 )l and m_(r) = m. 

To solve the system of equations (2.2) for the Green functions, we assume 
they are of the form G SlS2 = G® lS2 +Gl lS2 + G 2 lS2 where G\ lS2 are the solutions 
previously found for one impurity located at Ri . Without lost of generality we 
choose R 2 > R\. The space is now divided into three regions. The first region 
is n < the second region is R\ < r\ < R2 and the third one is R2 < r\. 
The differential equations satisfied by G± + are 

m[l + aiJ(r 1 -Ri) + a 2 <y(ri - R 2 )]G^_ + (r 1 ,r 2 ) 
+ma 2 S(r 1 - R 2 ) [G^ + ( ri ,r 2 ) + G^+(ri, r 2 )] 

+2d Zl G 2 _ + (r U T 2 ) = C 
2d- Zl G 2 ++ (r 1 ,r 2 )+mG 2 _ + (r 1 ,r 2 ) = C 



(3.17a) 
(3.17b) 



with general solutions of the form 



G^ + (r l5 r 2 ) 



Gi + (n,r 2 ) = { 



E ieZ e^A (3) ^( 



AG2 

E, 



1 ; 

J10! r(i) 



[mri j 



(2) 



n < i?i 

Ri < n < R 2 

R2 < n 

n < R2 
i?i < n < i? 2 
i? 2 < ri 

(3.18) 



The delta functions in Eq. (3.17a) impose that G 2 ++ is continuous at ri = Ri and 
ri = R 2 but G 2 _ + is discontinuous. Imposing these four boundary conditions, 
we find after some algebra 



(2) 



B, 



4-7r2?7 



m 2 a 2 

27T7? 

m 2 a 2 

27T77 

= B, (2) - 



Ii{mR 2 )e- ae ^ [l + am ]£>i + A 



(2) 



^(mi^e"^ [l + aino]Di, 



m 4 aia 2 -K'o(w|Ri — R2 

47T 2 ?7 



Ki{mR{)e- lWH ^ Di. (3.19) 
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with 



m m 3 aiir (m|Ri - r 2 |)if (m|Ri - R 2 |) 
L>\ = -— K Q (m\r 2 - R 2 ) H 2 r n , 

(3.20) 

4 

777 G^CK 2 2 

1 + aino + a 2 7to + otiot2n,Q — — -[^o(m|Ri-R 2 |)] . (3.21) 



The functions G\ + are therefore 



^ + (r!,r 2 ) 



G 2 _+(r 1; r 2 ) - 



m 2 a 2 
2ttt] 

m 2 a\ 
2ir 



m 2 a 2 

27T77 

m 2 ai 
2tt 



[l + ainojifoMri - R2I) 



#o(m|ri - Ri|)if (m|Ri - R 2 |) 



(3.22a) 



[l + aino]A-i(m|ri -R 2 |)e iei ^ 



A"i(m|ri - RiQe^iifoMRi - R 2 | 



.(3.22b) 



Following some similar steps, we find the other two Green functions 



G 2 __( ri ,r 2 ) = - 



m 2 a 2 

2-KTj 



[1 + am ] ifiHn - R 2 |)e iei «2 



m 2 ai 
2tt 



#i(m|ri -RiDe^itfoMRi -R 2 |) 



G+_(ri,r 2 ) 



m 2 a 2 

2-KTj 



Do 



[l + ain ]K (m\r-L - R 2 | 



m 2 ai 
2tt 



K (m| ri - RiDifo^IRi - R 2 |) 



(3.23a) 



with 



771 



D 2 = -—K 1 (m\R 2 -v 2 \)e- w ™2 



2ir 



+ 



m 3 a 1 K 1 (m\R 1 - r 2 |)i^ (m|R 2 - Ri|)e- 4e2 «i 



4tt 2 [1 



ain a \ 



(3.24) 



Replacing these Green functions in Eq. (2.3a) gives the density profiles. 
Figure 8 shows the charge density profile. It can be seen that the charge density 
decreases around the adsorbing points like it was found in the problem with 
one impurity. The density depends on the adhesivity as well as of the distance 
between the impurities. An important effect arises when the adsorbing points 
are close enough which indicates the existence of an indirect interaction between 
the impurities (see section 3.2.3). 
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Figure 8: Charge density profile for two "positive" impurities. The adhesivities 
are m^ai = 0.2 and m 2 a 2 — 0.5. The positions of the impurities are on the 
x axis at mR\ = 1, mR\ = 1.4 and mi? 2 = 2, mi? 2 = 1-6, respectively from 
the darkest to the lightest. All the angular differences are zero. The cutoff is 
ma = 0.01. 



3.1.3 One positive and one negative impurity 

Now we consider the case when the second impurity only adsorbs the negative 
particles. The fugacities for this configuration are m+(r) = m[l + ai<5(r — Ri)] 
and m_(r) = m[l + a 2 S(r — R 2 )] . Assuming a similar solution as in Eq. (3.18), 
and after establishing the appropriate boundary conditions for this case, we 
solve the linear system (2.2), to find 



G++(ri,r 2 ) 



G ( [ + (ri,r 2 )+GV + (r 1 ,r 2 



+ 



'++ 
m 2 a 2 
2nf] 

m 2 ai 
2tt 



D, 



[1 + amo^iMr! - R 2 |)e-^«2 



#oH r i - RiDe-^i^JTiHRi - R 2 | 



G_ + (r 1; r 2 ) = G°_ + (r 1 ,r 2 ) + GL+(r 1 ,r 2 ) 



+ 



m 2 a 2 
2nf] 

m 2 ai 
2tt 



[1 + amoje^JfoHri - R 2 |)e- 4ei «2 



(3.25) 
(3.26) 



^i(m|ri - RiDe^i-^i^^^mlRi - R 2 | 
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with 



D 3 = —e i6 ^K l {m\R 2 -v 2 \) 
Ztt 

m 3 a 1 K (m\'R 1 - T 2 \)e ieR i R 'K 1 (m\R2 - R^) 
4tt 2 [l + axTio] 

4 

2 m a\a 2 r 

?7 = 1 + ain + a 2 n Q + a^n^ -\ — - — |Ai(to|Ri 

and O^Ri the polar angle of the vector Ri — R2. 

In order to find G — , we note that if we interchange Ri and R2 and their 
adhesivities in the problem for G++, we would have a negative impurity at 
ri = Ri, and a positive one at ri = R2. Due to this symmetry argument, the 

function G for one positive and one negative impurity at ri = Ri and ri = 

R 2 , is the function G ++ with Ri and R 2 and their adhesivities interchanged. 
Then, 

G — (ri,r 2 ;Ri,ai,R2,a2) = G++(ri, r 2 ; R 2 , a 2 , Ri, ai). (3.29) 

with G ++ given in Eq. (3.25). 

Figure 9 shows the charge density profile for this situation. 



(3.27) 
-R 2 |)] 2 .(3.28) 



3.1.4 General structure of the Green functions for two impurities 

The results of the two previous sections can be put in more general and compact 
form. Let Si = ±1 be the sign of the particles the impurity at Ri attracts and 
S 2 = ±1 the sign of the particles the second impurity at R2 attracts. Let 
us denote by G SlS2 = G° lS2 + G\ lS2 , the Green functions for the case of only 
one impurity located at Ri (Eqs. (3.10) and (3.11)). Notice that in general 
Eqs. (3.10) and (3.11) can be written as 

G slS2 (r 1; r 2 ) = G° iS2 (r 1; r 2 ) - - — G° s s (r 1; R^Gg (Ri, r 2 ) 

(3.30) 

in terms of the Green functions G° iS2 and the density of the unperturbed 
system. 

The solutions for the Green functions found in the two previous sections for 
the case of two impurities can be written, in term of the Green function for one 
impurity, as 

G SlS2 (n,r 2 ) = G SlS2 (n,r 2 ) - ^G Sl s 2 (ri, R 2 )Gs 2S2 (R 2 , r 2 ) (3.31) 

where -q was defined in Eqs. (3.21) and (3.28) (denoted fj in this last equation). 
In general 

77 = l + ain Si (Ri)+a 2 n S2 (R 2 )+aia 2 n s ' i ( s ) 2 (Ri,R 2 ) (3.32) 
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Figure 9: Charge density profile for one "positive" impurity and one "negative" 
impurity. The adhesivities are m 2 a\ = 0.2 (positive impurity) and rr?ai = 0.5 
(negative impurity) . The positions of the impurities are on the x axis at mR\ = 
1, mRi — 1.4 (positive impurity) and mi?2 = 2, to_R 2 = 1-6 (negative impurity), 
respectively from the darkest to the lightest. All the angular differences are zero. 
The cutoff is ma = 0.01. 
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where n^^g (Ri, R 2 ) is the (non-truncated) pair correlation function for the un- 
perturbed system. Notice the similarity in the structure of the Green functions 
in Eqs. (3.30) and (3.31). 

Replacing Eq. (3.30) into Eq. (3.31) gives the expression 



max [1 + a 2 n° S2 (R 2 )]G° lSl (n , Ri)G° SlS2 (R 1; r 2 ) 



G SlS2 (ri,r 2 ) = G° lS2 (ri,r 2 ) 
1 
V 

+ma 2 [l + ai n° Sl (Ri)]G° iS2 (n, R 2 )G° 2S2 (R 2 , r 2 ) 
-m 2 ai a 2 [G° SiSi (n, Ri)G° SiS2 (R 1; R 2 )G° 252 (R 2 , r 2 ) 

+G2 lS2 (r 1 ,R 2 )G° S2Sl (R 2 ,R 1 )G° lS2 (R 1 ,r 2 )] 

(3.33) 

in which the exchange symmetry between (Ri,Si) and (R 2 ,S' 2 ) is manifest. 



3.2 General results 

The grand potential of the system with impurities can be found by solving the 
eigenvalue problem (2.6) for each particular situation. We will proceed this 
way to compute the grand potential for the Coulomb gas in the presence of a 
continuous adsorbing line in subsections 4.2.1 and 4.2.2. Although this method 
can also be used to compute the grand potential with point impurities [13], for 
this latter case, we shall proceed differently, using the general theory developed 
in Ref. [8] for a one-component fluid with adsorbing sites. 

In the following subsections, we obtain expressions for the grand potential 
and the one- and two-body densities of the plasma with impurities in terms of 
those same quantities for the plasma without impurities (the unperturbed sys- 
tem). These results are very general: they are valid for any value of the coupling 
T, they are even valid for any other type of interaction potential between the 
particles of the system, not only for the Coulomb system considered here. In 
the present case, since at T = 2 the grand potential and density functions for 
the unperturbed system are known, one can obtain explicit expressions for the 
thermodynamics and correlations of the plasma with impurities. 



3.2.1 Grand partition function with an arbitrary number of impuri- 
ties 

Suppose there are M + impurities attracting positive particles, located at R+ , . . . , R ^ 
with adhesivities a^, . . . , a\ I+1 and M_ impurities attracting negative particles, 
located at R^, . . . , R^ f _ , with adhesivities af, . . . , ot~ M _ . The general calcula- 
tions done in Ref. [8] for a one-component plasma, can directly be adapted to 
the present two-component case, to find that the grand partition function 5 can 
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be written as 



oo oo M+ M_ TT S + „ T\ s 

S °E E E E ^' n '- n "-'' 



s+!s_! 



s + =0 s_ =0 fe+ , =0fe~ ,...,fe7_ =0 



where ^o is the grand partition function of the unperturbed system and n°^^ s s ? 
the (s + + ,s_ )-body density for s + positive particles and s_ negative particles 
of the unperturbed system. Notice that, since the correlation function vanishes 
if two of its arguments are equal, the above expression is not an infinite sum, 
it contains at most terms involving the (M + + M_ )-body correlation function 
and lower degree correlations. 

From expression (3.34) we can obtain the density and fc-body correlation 
functions performing successive functional derivations of the grand partition 
function with respect to the fugacity. 

3.2.2 One impurity 

In the case of a single impurity located at Ri, with adhesivity cti and attracting 
particles of sign Si, Eq. (3.34) simply reduces to 

S= [l + am Sl (Ri)]So. (3.35) 

Then, the grand potential £1 = — fc^TlnS can be written as Q = Qq + f2 cxc (ai) 
with ilo the grand potential of the unperturbed system and an excess grand 
potential 

/3ft e xc(ai) = -ln[l + ain Sl (Ri)] . (3.36) 

As explained in Ref . [8] , the adhesivity a can be interpreted also as a fugac- 
ity for the adsorbed particles. Thus, one can compute the average number of 
adsorbed particles from the relation 

Ng 1 = _ Ql g^£M = . (3.37) 

Sl da! l + am^Ri) v ' 

We recover the result (3.14) obtained from a direct calculation at V = 2. 

From Eq. (3.35) we can also rederive the results from the previous section 
for the density and correlation functions, by performing functional derivatives 
of S with respect to the fugacity. Remembering that 



(HnSo 

m s (r) 



Sm s (r) 

and 



= n° s (r) (3.38) 

i(r)= m 



m S2 (r 2 ) 



5m S2 (r 2 ) 



= <(f>i,r 2 ), (3.39) 

m(r)= m 
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deriving (3.35) with respect to m s (r), we obtain 

0,(2)T, „ x 



n Sl (r) = [l + ai<y(r-Ri)] 



n° Si (r) 



and 



n-Si(r) = n°. Sl (r) 



l+ai4( R i) 

0.(2)T/ -r, \ 



(3.40a) 



(3.40b) 



l + am° Si (Ri) 

When r = 2, we recover the results (3.12) from section 3.1.1. Deriving Eqs. (3.40) 
once again with respect to the fugacity, we obtain the truncated two-body cor- 
relation functions 

[1 + oi<y(ri - Ri)] [1 + ai<5(r 2 - R^] ng^*(ri,r 2 ) 

(3.41a) 

[1 + ai S(r 2 - Ri)] (n , r 2 ) (3.41b) 

(3.41c) 



( 2 ) T ( \ 

n^ Sl (ri,r 2 ) 



( 2 ) T I \ 

with 



,(2)1 



^ li _S 1 (ri,r 2 ) 



n 



(2)T* 



0,(2)T, 



(ri,r 2 ) 



(ri,r 2 ) 



+ 



0,(2)T 

0,(3)T , D \ 



ain^(Ri) 



(3.42) 



( 3)T 

where n s [ S 2s 3 is the truncated three-body density function of the unperturbed 
system. From the Green functions, Eq. (3.30), obtained in the previous section 
and the general relations (2.3), it can be shown that the correlation functions 
obtained in section 3.1.1 are indeed those given by Eqs. (3.41) and (3.42) when 

r = 2. 

3.2.3 Two impurities 

Now, let us consider the case of two impurities, located at Ri and R 2 with 
adhesivities cti and a 2 and attracting particles of sign Si and 52, respectively. 
Eq. (3.34) gives 

S = S [l + a 1 n^ i (R 1 ) + a 2 n^ 2 (R 2 ) + aia 2 n°' i ( |(R 1 ,R 2 )] = E oV (3.43) 

with r/ defined in Eq. (3.32). The grand potential now reads 

0Q = /3Q -In 77. (3.44) 

Notice that in the problems for one and two impurities the excess grand potential 
is the logarithm of the denominator of the Green functions G], lS2 and Gg lS2 
respectively. It is convenient to express the grand potential as 



fl = Q + fiexc(ai) + ^0x0(^2) + ^5iS 2 ( R l7 «i; R 2, OI2) 



(3.45) 
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Figure 10: Effective interaction between two positive impurities as a function of 
their distance. The adhesivities are m 2 a\ = m 2 a: 2 = 0.2, 1 and 2 from darkest 
to lightest. The cutoff is ma — 0.01. 



with 



/^SiS 2 (Ri,ai;R-2,a;2) = - In 
Explicitly, for T = 2, 

/3fi ±± (Ri,ai;R 2 ,a 2 ) = - In 

and 

/30 ±=F (Ri,Q!i;R2,a2) = - In 



aia 2 n s ^ s ' 2 (Ri,R 2 ) 



1 + 



[l + am^(Ri)] [l + a 2 n° 2 (R 2 )] 



m 4 aia 2 [i4T (m|Ri -R2I)]' 
47r 2 [l + aino] [l + a 2 n ] 



1 + 



n 4 o:ia 2 [ifi(m|Ri — R 2 |) 
47T 2 [l + aino] [l + a 2 n ] 



(3.46) 



(3.47) 



(3.48) 



A work n exc (a>i) is required to introduce a single impurity into the plasma. 
To introduce two impurities, we require, in addition to O exc (ai) + r2 0XC (ai 2 ), 
an additional work 05 1 s 2 (Ri, ai; R 2 , a 2 ). We can define the effective chemi- 
cal potential of an impurity with adhesivity a as n{a) = fl cxc (a). The term 
Os 1 S 2 (Ri,ai;R2, a 2 ) represents an effective interaction between the impurities. 
This contribution is not a direct interaction between the adsorbing particles but 
a consequence of the inhomogeneous charge distribution that they create. 
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Figures 10 and 11 show the effective interaction between two adsorbing par- 
ticles. It can be seen that it is repulsive if the impurities have the same sign and 
attractive if their sign is opposite. As expected, the interaction is important 
when the impurities are close enough but is negligible for large distances. It 
depends on the adhesivity of the two impurities, indeed, it is larger for large 
adhesivities. 

From Eq. (3.43), we can obtain the density profiles by performing a func- 
tional derivative with respect to the fugacity. We obtain 



n s (r) 



'[l + aiJ(r-Ri)+a 2 J(r-R2)]n;(r) if s = Si = S 2 

[l + aiJ(r-Ri)]n;(r) if s = Si = -S 2 

[l + a 2 J(r-Ra)K(r) if s = S 2 = -Si 

IX (r) if s = -Si = -S 2 



(3.49) 



with 



ain^foRi) [1 + a 2 n° S2 (R 2 ) 



+ a 2 n s jl )T (r, R 2 ) [l + am° Sl (Ri)] + aia 2 u°^ (r, Ri, R 2 



0,(3)T / 



(3.50) 

With the Green functions (3.33) for two impurities found in sections 3.1.2, 3.1.3 
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and 3.1.4 and the general relation (2.3), it is direct to verify that the density is 
indeed given by Eqs. (3.49) and (3.50) when F = 2. 

We can compute the average number of adsorbed particles by each impu- 
rity Ng* , for i = 1,2. This can be obtained either from the relation Ng* = 
ajng.(Rj), or from 

N ^ = - ai -^T- (3 ' 51) 



We obtain 

N st = 



<Xjn Si ( R + a ^ n °s[ 2 S 2 ( R i ; R 2) ^ 

1 + ain^^Ri) + a 2 n S2 (R 2 ) + ai^^^Ri, R 2 ) 



4 The plasma with adsorbing lines 

4.1 Density and correlations 
4.1.1 An infinite line 

As a simple model for electrodes with adsorbing sites, we consider now that a 
line of impurities is introduced into the plasma. We shall consider the case where 
the impurities are very close to each other, so there is a continuous adsorbing 
line in the plasma, located at x = x . Since the system is translational invariant 
in the direction of the line (the y direction), it is better to work in Cartesian 
coordinates and confine the plasma in a rectangular box. The fugacity can be 
modeled by 

m ±( r ) = m±(x) — m[l + a±5(x — x )] ■ (4.1) 

The adhesivity a± has length dimensions, unlike the adhesivity a± for point 
impurities which has area dimensions. Actually, a± is an effective adhesivity 
due to the distribution of the adsorbing particles on the line. We will assume 
that the line is composed of "positive" impurities, so a + — u\ and a_ = 0. 

Due to the translational invariance in the y direction, this problem is easier 
to solve via Fourier transform. For this purpose we define 

G aiaa (ri,r 2 ) = -L- f G SlS2 (x 1 ,x 2 ,l)e il ^-^dl. (4.2) 



By replacing Eq. (4.2) and the given fugacities into Eqs. (2.2), we obtain the 
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system 

m[l + ai6(xi - x )]G ++ (x 1 , x 2 , 1) 

+ (^-+l S jG-+{x 1 ,x 2 ,l) = 6(xi-x 2 ) (4.3a) 

(-^-f)G ++ (x 1 ,x 2 ,l)+mG-+(x 1 ,x 2 ,l) = (4.3b) 
m[l + ai<5(xi - x )]G+_(xi, x 2 , Z) 

+ (^-+z)g__(xi,x 2 ,Z) = (4.3c) 

(dx ')G+-(ii|i2 ) i)+ m G ! — (xi,x 2 ,l) = 5(xi—x 2 ). (4.3d) 

We assume solutions of the form G S1S2 = G° iS2 + Gl iS2 where G° iS2 are the 
bulk solutions for a.\ — 0. It is easy to show that these solutions are given by [1] 



G° ±± (x u x 2 J) = ^e- k \*>-*>\ (4.4a) 
G°_ + (x u x 2 ,l) = ±(--jL+f)G ++ (x u x 2 ,l) (4.4b) 
G° + _{ Xl ,x 2 ,l) = -±(-jL+i)& i __(x 1 ,X2,l) (4.4c) 



with k = V™ 2 + I 2 - 

The solutions to Eqs. (4.3) are 



G ++ (xi,x 2 ,l) 



t e- fc i— i - 2fe(2 r + ^ ai) e " fc|23;o " xl " 3;21 ' 2:0 < 

or x > xi,a; 2 (4.5) 

-fe|xi— X2| 



=-e |:El 121 x between x\ and x 2 



v 2k+m 2 ai 

!m -k\x!-x 2 \ , ma^fc+Q 2 fc(2x -si-x 2 ) ^ 
2fc+m^ 1 e " fc|Xl " X2 ' > between an and x 2 (4.6) 

m p -fc|xi-x 2 | 1 ™fii(fc-Q 2 fc( 2xo _ Xl _ X2 ) 
2fe e + 2fe(2fc+m 2 ai) e ' X ^Xi,X 2 

The number density profiles n+, n_, and charge density profile p (in units 
of e) are obtained using Eqs. (2.3). We find 

n + (x) = [1 + a\5(x — xo)] n* + {x) (4.7a) 
n-(x) = n*_(x) (4-7b) 
p(x) = n+(x) — n-(x) = a\5(x — xa)n* + (xo) + p*(x) (4.7c) 
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Figure 12: Positive density for a line of "positive" impurities at mx — 1. The 
adhesivities are mdi — 0.2, 1 and 3 from the darkest to the lightest. The Dirac 
distribution at mx — 1 is not plotted. The cutoff is ma — 0.01. 



with 

n*Ax) = n°, - m 3 ax / , — , (4.8a) 

+v ' + Jo 47rVFTT(VFTT + a) v ; 

r°° dt C1 4- It 2 ) p -2Vt 2 +lm\x-x \ 

n*(x) = n° +m 3 a 1 — (1 t_L 6 = = (4.8b) 

p {x) = - m a ^L ^ vFTTTa (48c) 

with a = ma.\j1. 

Figures 12,13 and 14 show the positive and negative density and the charge 
density profile for this configuration. The curves obtained are similar to the 
ones obtained for point impurities. They differ in the magnitude as well as in 
the peaks near the adsorbing particles. The magnitude of the density at the 
adsorbing line is reduced as a consequence of the impurities distribution. 

The linear charge density adsorbed by the line, a, can be computed by the 
relation 

o = a in * + {xn) (4.9) 
or by integrating the opposite of the nonadsorbed charge density 

/•oo 

<j = - I p*(x)dx. (4.10) 
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Figure 13: Negative density for a line of "positive" impurities at mx = 1. The 
adhesivities are ma.\ — 0.2, 1 and 3 from the darkest to the lightest. The cutoff 
is ma = 0.01. 
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Figure 14: Charge density profile for a line of "positive" impurities at mx = 1. 
The adhesivities are ma\ — 0.2, 1 and 3 from the darkest to the lightest. The 
cutoff is ma = 0.01. 
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Either way, the result is the same, as expected, 



m 2 a x ^ ^ ,,, '., , . (4.Ha) 



dt 



n m a 2 \Ja 2 — 1 
= ainrL arctanh (4.11b) 

where we introduced an ultraviolet cutoff t max = e~ c /(ma) to obtain a finite 
result. This choice of the cutoff is done to ensure the correct result for bulk 
density [1]. 

Contrary to the case of a point adsorbing impurity, where the adsorbed 
charge reached a maximun value of 1 when a\ — > oo, the adsorbed line charge 
density is not bounded when a\ — > oo. Indeed, for large <5i, 



n ma 

OL\Vb\ 

7T 



ln(2a) + ±- 2 



In(2a) - \ 



+ 0(a- 4 ) 



(4.12) 



It is instructing to compute the electric field E(x) and electric potential 
<p(x) created by the adsorbed line charge density and its polarization cloud. 
Integrating Poisson equation we obtain, for x > and choosing x = 0, 

/■oo -2\/t' 2 + lmx ji 

E(x) = ema T^f—- (4-13) 

Jo Vt + 1 + a 

and 

4>(x) = ea / (4.14) 

Notice that the electric field at the position of the impurity is related to the ad- 
sorbed linear charge density by E(0 + ) = ncr, and as usual we have the expected 
discontinuity E(0 + ) — E(0~) = 2tt<j. The potential drop from the adsorbed line 
to infinity is 



(f)(0) - 4>(+oo) = ea [ 
Jo 



dt ea \Ja 2 — 1 
: arctanh 



t 2 + 1 + aVt 2 ^ Va 2 - 1 " 

^ (4-15) 

The correlation functions can be obtained from the Green functions Eq. (4.5) 
and (4.6), using the relation (2.3b). Let us simply note that the correlation 
functions have an exponential decay both in the x direction (transverse to the 
adsorbing line) and the y direction (along the adsorbing line) , a behavior similar 
to the case without adsorbing line. This is different from the algebraic decay 
that the correlation functions show parallel along a hard wall [14, 15, 16]. 

4.1.2 Two infinite lines 

For two parallel infinite lines it is also possible to find analytic expressions for the 
Green functions. For two lines located at x = X\ and x — X 2 , both attracting 
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the positive particles, the fugacities are 

m+(x) = m[l + a 1 S(x-X 1 ) + a 2 S(x-X 2 )} (4.16) 
m-(x) = m (4.17) 

For one line, located at x = X\, attracting the positive particles, and a second 
line, located at x = X 2 , attracting the negative particles, the fugacities are 

m+(x) = m(l + an6(x - X^) (4.18) 
m-(x) = m(l + a 2 S(x - X 2 )) (4.19) 

Without loss of generality, we suppose that X\ < X 2 . The resolution of the 
linear system (2.2) satisfied by the Green functions follows similar steps as for 
the previous case of one line. The explicit expressions for the Green functions 
depend on the relative positions of their arguments with respect to the lines. 
The complete expressions can be found in Rcf. [13]. Let us focus here in the 
expressions obtained for the density profiles. 

For two lines attracting positive particles, we find the positive, negative and 
charge densities 

n+(x) = [l + a 1 6(x-X 1 )+a 2 S(x-X 2 )]n* + (x) (4.20a) 
n-(x) = n*_{x) (4.20b) 
p(x) = [a 1 S(x-X 1 ) + a 2 S(x-X 2 )]n* + (x)+p*(x) (4.20c) 

with 

f 00 m 4 e- 2kx a 1 {2k-m 2 a 2 )e 2kXl +a 2 (2k + m 2 a 1 )e 2kX2 „ ,,, „, , 

n,(x) = n, — I — - — t-, s — ; 3 ttt^ — ^rr dl 4.21a) 

+ K ' + i^oo Airk (2k + m 2 a 1 )(2k + m 2 a 2 )~m 4 aia 2 e 2k ^- x ^ V ' 

n f°° m 2 (k + l) 2 e- 2kx aA2k - m 2 a 2 )e 2feXl + a 2 (2k + m 2 aAe 2kx * „ 

~ V ; " 4irk (2k + m 2 ai)(2k + m 2 a 2 ) -m 4 aia 2 e 2k< - x i- x ^ 

(4.21b) 

f + °° m 2 ke- 2kx ai(2fc - m 2 a 2 )e 2kXl + a 2 (2k + m 2 ai)e 2fcX2 
9 W ~ 'Loo 2^ (2k + m 2 a 1 )(2k + m 2 a 2 ) - m 4 a x a 2 e 2k ^- x ^ ( ' C) 



for X\ < X 2 < x, with k = yjm 2 + I 2 . Between the lines, we find 

n* + (x) = n + 



I 



m 4 on(2k + m 2 a 2 )e 2k(Xl - x ) + a 2 (2k + m 2 a 1 )e 2k ^- X2 ^ - 2m 2 a 1 a 2 e 2k( - Xl - x ^ 



00 



47rA: (2k + m 2 ax)(2k + m 2 a 2 ) - m 4 a 1 a 2 e 2k( ~ Xl - x ^ 

(4.22a) 

n*_(x) = n°_ 

m 2 ai(fc + l) 2 (2k + m 2 a 2 )e 2k ^- x ^ + a 2 (k - l) 2 (2k + m 2 ai)e 2fc ( x - X2 > + 2m 4 a 1 a 2 e 2k{ - x ^ x ^ 
4irk (2k + m 2 a 1 )(2k + m 2 a 2 )~m 4 aia 2 e 2k< > Xl - x ^ 

(4.22b) 

f + °° m 2 k 6n(2k + m 2 a 2 )e 2k ^-^ + a 2 (2k + m 2 ^ 2 ^-^ 
P{X) ~ ~ 7-oo "2^ (2k + m 2 a 1 )(2fc + m 2 a 2 ) - m 4 aia 2 e 2fe ^-^) ( ' 
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for X\ < x < X 2 . It is straightforward to verify that these expressions reduce 
to case of one adsorbing line when X\ — > — oo, as it should. Also, one can verify 
that for X\ = X 2 , we recover the density profiles for one line with adhesivity 
a\ + a 2 . Figure 15 shows the charge density profile in the presence of two 
positive lines. 

The linear density of adsorbed particles by each line, o\ and a 2 , are obtained 
from u\ — din* + {Xi) and a similar expression for a 2 . Explicitly, 

m 4 ai f + °° [a l (2k + m 2 a 2 )+a 2 (2k-m 2 a 1 )e 2k( - x ^- x ^] dl 
ax - a x n + — J ^ - - m2 -^ 2k ~ m 2 a ^ _ ^a^e 2 ^-*^] 

+°° (2k + m 2 a 2 )a x - m 2 ^^ 2 ^^ 

(2k + m 2 a 1 ){2k + m 2 a 2 )-m i a 1 a 2 e 2k ^- x ^ ^' ' 

and a similar expression for a 2 interchanging ck\ and a 2 . In the second line of 
Eq. (4.23) we used the formal expression 

for the bulk density (actually, the integral should be cutoff to a Z ma x = e c /a 
to obtain finite results). The total adsorbed charge is found to be 



to* r 
^r" J_ 



CTl+fJ2_ J-oo (2k + m 2 a 1 )(2k + m 2 a 2 )-m i a 1 a 2 e 2k ( x ^ x ^ ( ' } 
It is straightforward to verify that 

r+oo 

0-1+0-2 = - / p*(x)dx (4.26) 



as it should be, since the system is globally neutral. 

In the case where line at X 2 adsorbs negative particles while the line at X\ 
adsorbs the positive ones, we find the following density profiles 

n+(x) = [1 + ondix - X-C)]n* + (x) (4.27a) 
n-(x) = [1 + a 2 5(x- X 2 )]n*_(x) (4.27b) 
p(x) = a 1 S(x-X 1 )n* + (X 1 )-a 2 S(x-X 2 )n*_(X 2 )+p*(x) (4.27c) 



with 



. . . 1 1 / 772 e ni uni^ — /it U2;c ~ — ix 2 \n> — i I lift t "t Ul ic 

ji I x) = n — / - - - - - - dl 

' " ! 4irk (2k + m 2 a 1 )(2k + m 2 a 2 )+m 2 a 1 a 2 {k + l) 2 e 2k ( x ^- x ^ 

(4.28a) 
(4.28b) 

» m 2 e -2kx ^ k + t ^ 2k _ m 2^ e 2kx 1 + ^ _ k ^2k + m'aQe 2 ** 2 
2^ (2fc + m 2 ai)(2fc + m 2 a 2 ) + m 2 (k + l) 2 a x a 2 e 2k ^- x ^ 

(4.28c) 



-2kx m 2^2k - m 2 a 2 )e 2kXl - a 2 (k - l) 2 (2k + m 2 ai)e 2kX2 



„ . m 2 e" 2fex ai(fc + 2 (m 2 a 2 - 2fc)e 2fcXl +m 2 a 2 (2fc + m 2 a 1 )e 2 ' cX2 
~ n -~ I > 4 ^ fc (2fc + m 2 ai)(2fc + 777 2 a 2 )+m 2 aia 2 (fc + 2 e 2fc ( Xl - X2 ) 
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m 2 [n+ -n_ ] 

2r 



1 - 




Figure 15: Charge density profile as a function of the distance. The adhesivitics 
are ma>i = 0.2 and mo.2 =0.3. The position of the lines are mX\ = 1, mX\ = 
1.3 and 111X2 — 2, 111X2 = 1.7 from the darkest to the lightest. The cutoff is 
ma = 0.01. 



for Xi < X 2 < x. For x between the lines, X\ < x < X 2 , we find 



n* + (x) = n° + 

r°° m 2 m 2 a 1 (-2fc - m 2 a 2 )e 2k ^-^ + a 2 (k + l) 2 {2k + m^Qe 2 ^^) - 2m 2 a 1 a 2 (k + QV^-X^ 
+ J_ co 4^k (2fc + m 2 a 1 )(2k + m 2 a 2 ) + m 2 (k + l) 2 a x a 2 e 2k ^- x ^) 

(4.29a) 

n*_(x) = n°_ 

r°° m 2 a x (k + Q 2 (2fc + m 2 a 2 )e 2fc(Xl ~ a:) - m 2 a 2 (2fc + mV)e^^ ~ 2m 2 aia 2 (fc + ^V^-^) 
+ i_ 00 i^fc (2fc + m 2 a 1 )(2k + m 2 a 2 ) + m 2 (k + Ifoa^e 2 ^-^ 

(4.29b) 

f +00 m 2 (k + l) q 1 (2fc + m 2 Q 2 )e 2fc(Xl ~ :E) -a2(2fc + m 2 ai)e 2fc < a: ^ X2 ) 
1 (X) ~ J.^ 2^ (2fc + m 2 a!)(2fc + m 2 a 2 ) + m 2 (k + l^a^e 2 ^-^ ' ^ ' 

The linear charge density adsorbed by each line is now given by o\ = 
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ain* + {Xi) and a 2 = 



-a 2 n*_(X 2 ), 



2 r+°° 



aim 2 r 
2tt J_ 



(72 = 



[2k + m 2 a 2 + {k + l) 2 a 2 e 2k{x ^ x ^} dl 
(2k + m 2 ai)(2fc + m 2 a 2 ) + m 2 (k + l) 2 ai a 2 e 2k< - x ^- x ^ 

(4.30a) 

a 2 m 2 f +oc [2k + m 2 a 1 + (k + l) 2 a 1 e 2k< ^ Xl ~ x ^] dl 

2tt J_ 00 (2k + m 2 a 1 )(2k + m 2 a 2 ) + m 2 (k + l) 2 a 1 a 2 e 2k ( Xl ~ x ^ 

(4.30b) 



The total linear charge density adsorbed by the lines is formally given by 
^ 2 r +co kdl 



m f 
= (ai— a 2 ) — / 

7T J-< 



[2k + m 2 a 1 )(2k + m 2 a 2 ) +m 2 (k + l) 2 a x a 2 e 2k ^- x ^ 

(4.31) 

Once again, it is straightforward to verify that 



<J\ + <J 2 



p*(x)dx (4.32) 



as imposed by the global neutrality of the system. 

It is interesting to notice that in the two cases considered in this section, 
the density profile of particles that are adsorbed by one line is a continuous 
function as it crosses the line, while the density profile of the particles that are 
not adsorbed by this line is discontinuous when crossing such line. However, the 
discontinuity jump in this density is exponentially small for large separations 
between the two lines, it is of order e - 2m \ x i- x *\ . 

4.2 Grand potential 

In the following subsections, we compute the grand potential of the plasma with 
one or two adsorbing lines at T = 2, by solving the eigenvalue problem (2.6). 

4.2.1 An infinite line 

For this case we confine the plasma into a rectangular box of area 2L X x L y 
and we work in Cartesian coordinates. The boundaries and the adsorbing line, 
located at x = 0, divide the space in four regions x<— L x , — L x < x < 0, 
< x < L x and x > L x which will be labeled by (1), (2), (3) and (4). 
The eigenvalue system (2.6) takes the form 

m[l + ai8(x)]g{r) = 2Xd z f{r), (4.33a) 
m/(r) = 2Xd s g{r). (4.33b) 

for — L x < x < L x . Outside that region, m + (x) = m-(x) = 0. Thus we 
conclude that for x £ [—L X ,L X ], d z f = and d s g = 0, that is / is an anti- 
analytic function of z and g is an analytic function of z. 
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The translation symmetry on the y axis allows us to assume g(r) of the form 
g(x,y) = g{x)e lly . The general solution for the four regions can be written as 



g {2) (x,y) = AWe k * x+ih - 
9 {4) (x,y) 



B {2) e -{k xX - 



^(3) e k x x+ily _|_ £j(3) e -(k x x-ly) 



and 



f (1A \x,y) = C^e- lx+ily 

/ (2 ' 3) (*,</) = 



d 



dy 



,(2,3) 



{x,y). 



with k. 



-ft 



+ ■ 



(4.34a) 
(4.34b) 
(4.34c) 
(4.34d) 

(4.35a) 
(4.35b) 

(i) = 



Since the eigcnfunctions must vanish at x — ±00, we conclude that g 
/( 4 ) = for I < and /W = g^ = if I > 0. The boundary conditions 
demand that both g(x, y) and /(a;, y) must be continuous at x = ±L X . At x = 0, 
g(x,y) is continuous and f(x,y) is discontinuous due to the Dirac distribution 
in Eq. (4.33a): 



/(0+ tf )-/(0-,y) 



ff(0,y). 



(4.36) 



Assuming I > 0, these boundary conditions can be expressed by the linear 
system 



/ (k x 


— l)e kxLx 


-{k x + l)e k * L * 


\ 


( A 2 


\ 







e fc - L - 


e -k xLx 


B 2 




1 


1 -1 


-1 






\ m 2 a\ 


-X 2 (k x -l) 


m 2 ai + X 2 (k x + I) X 2 (k x - I) 


-X 2 (k x + l) J 


V B 3 


J 



(4.37) 

The functions g(x,y) and f{x,y) must not be zero, hence we demand that the 
determinant of the latter matrix must vanish. This condition leads us to the 
relation 



cosh(2k x L x ) 



1 + 



m 2 a\lX 



-2-, 



2k 2 



+ simi^fe^) 



I m 2 ai\ 
1 — 

k x 2k x 



m 2 ail\ 



= 0. 



(4.38) 

For each value of I, there are several possible solutions to Eq. (4.38) for the eigen- 
value A, which we will denote as {X^ n }n- For I < 0, we obtain (4.38) changing 
I by — I. To obtain the grand potential, we can recognize its relationship with 
a Weierstrass product [18] as follows. Let us define the analytic function, for 
I > 0, 



hi(z) 



1 + 



+ 



2 - ; 2 

m a\lz 
2(m 2 z 2 + I 2 ) 

T72 



) cosh(2 V / m 2 z 2 + l 2 L x ) + 



m 2 a\lz 2 



m^z- 



00 



21 



2{m 2 z 2 + I 2 ) 
—J smh(2^m 2 z 2 + l 2 L x ). (4.39) 



0. 
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Wehave/i ; (0) 



, /i;(0) = and hi(z) = hi(—z). By construction, the zeros 



of hi are precisely the inverse of the eigenvalues 1/ A; 
hi(z) as a Weierstrass infinite product is therefore 



The representation of 
(4.40) 



h l {z) = h l {Q)\{{l-zX hn ) 

n 

Then, we notice that the grand potential is simply 

/30 = -VlnTT(l + A i . Il ) = -^ / dl In [e-^tni-l)] (4.41) 

where we replaced the sum over I by an integral in the thermodynamic limit 
L y — > oo. In the thermodynamic limit, we also consider that L x — > co, and 
we approximate cosh(2fcL a .) <~ smh(2kL x ) <~ e 2feix /2, with fc = \Jl 2 + to 2 . The 
final result is 



/3ft = /3ft +/3fte X c(ai), 



with 



+00 



In 



TO 2 o:i 



2Vm 2 + P 



(4.42) 



(4.43) 



and fi the grand potential of the unperturbed system for this geometry, 

fl = -2L x L yPb + 2L yl + 0(e- 2mL * ) (4.44) 
with the bulk pressure pb and the surface tension near an impenetrable wall 7, 



a m 
PPb - 2^ 

Pi = 1 



In A 
ma 
1 " 

~ 2tt- 



TTi 

H0 + ^ 



(4.45) 
(4.46) 



Hence we can conclude that the effective chemical potential of a line of length 
L y and adhesivity a is fj,(a) = Clf xc (a). The integral (4.43) can be computed 
explicitly. Actually, one must introduce an ultraviolet cutoff Z max = e~ c /a to 
obtain a finite result. This choice of the cutoff is done to ensure that the result 
for bulk pressure (4.45) is the same as the one obtained by integrating the bulk 
densities (3.7) [1]. Neglecting terms that vanish when a — > 0, we have 



/JO exc (a) = - mL y 



7T 



a In ■ 



2e" 



UlCL 



— V^ 2 



— \J OL 2 — 1 7T 

1 arctanh ; — + a 

a 2 



(4.47) 



with a — ma/2. It is interesting to compare this result with the surface tension 
obtained in a similar problem where the plasma is confined in a strip with 
adsorbing boundaries [10]. In that problem the adsorbing line was located at a 
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hard wall boundary. In our notations, the contribution of both the hard wall 
and the adsorbing line was [10] 



/□(-^adsorbing wall UlLy 

P"cxc - 



-, 2e ~ c , - 1 - (™5) 2 , ,, 

ma In h 1 — n + ma H = ln( 1 + ma) 

ma a 



(4.48) 

We notice that £l exc (a)+L y j of our present problem is different from ri^ orblng wal1 . 
This is to be expected since each situation is different. In our present problem 
the adsorbing line is very far from the hard wall boundaries, while in Ref. [10], 
the adsorbing line is on the boundary itself. We notice however that the (di- 
vergent) dominant term, when the cutoff a vanish, for our problem is twice the 
one for the adsorbing boundary 

n cxc (a) ~ 2^ c sorbing wal1 ~ ^ ln(ma) . (4.49) 

The factor 2 between both expressions can be understood if we realize that in 
our present problem the plasma is on both sides of the adsorbing line, while in 
Ref. [10] it is only on one side, the other side is empty. 

From Eqs. (4.43) and (4.11a) or from Eqs. (4.47) and (4.11b), one can check 
that the relation between the adsorbed charge density a and the excess grand 
potential 

a = _ & ewiL y) 

oa 1 

is satisfied. 



4.2.2 Two infinite lines 

Let us now consider the case where there are two infinite parallel adsorbing lines 
which attract the same kind of particles, say positive. The eigenvalue problem 
which must be solved is 



m/(r) = 2Xd- z g(r) 
m[l + a 1 8(x-X 1 ) + a 2 (x~X 2 )]g(r) = 2A0 z /(r). 



(4.51a) 
(4.51b) 



We have placed the first line at x = X\ and the second one at x = X 2 - The 
method of solution is similar to the one for one line. Following similar steps as 
in the previous section, we find after some algebra [13] 



with 



0SI 



±± 



(3fl = f3n a + /3fi cxc (ai) + /30= X c(a 2 ) + /?^ ±± 
m i a 1 a 2 e- 2k \ Xl - x ^ 



— i f k 

T JO 



1 - 



(2k + m 2 a 1 ) {2k + m 2 a 2 ) 



dl 



(4.52) 



(4.53) 



and /3f2 e xc(a) given in (4.43). We recall that k — yjm 2 + I 2 . 
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Figure 16: Effective repulsion between two "positive" lines as a function of their 
distance 2xq = \X± — X 2 \. The adhcsivities are max = ma 2 = 0.2,0.5, 1 from 
the darkest to the lightest. The cutoff is ma = 0.01. 



When the lines attract different types of particles the grand potential is 

(3fl = f3Qo + /3f7 0XC (ai) + (3n cxc (a 2 ) + [3^ (4.54) 

with 

t r+°° 
f3Sl±T = _±v In 

The effective interaction between two lines (depending on their sign) is given 
by the expressions (4.53) and (4.55). 

Figures 16 and 17 show the effective interaction between two infinite lines. 
As expected, like in the problem for point impurities, there is a repulsion if the 
sign of the lines is equal and a repulsion if it is different. 

Notice that in both cases, s = ±, s' = ±, in the expression of the full 
excess grand potential f2 exc (ai) + ^6x0(0^2) + Q ss appears the logarithm of a 
term which is precisely the denominator in the Green functions and the density 
profiles Eqs. (4.21), (4.22), (4.28) and (4.29). Also, one can easily check that 
the adsorbed charge on each line, computed in Eqs. (4.23) and (4.30), can also 
be obtained from 

= -J^ a ^ (*W«i) + ^cxc(a 2 ) + ft 85 ') (4.56) 

for i = 1,2. 



m o.\o.i{k + I) e 



•2 -2fc|Ai-X 2 | 



(2k + m 2 ai) (2k + m 2 a 2 ) 



dl. 



(4.55) 
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Figure 17: Effective attraction between a positive and a negative line as a func- 
tion of their distance 2xo = \Xi — X2\. The values for the different parameters 
are the same as in figure 16. 



5 Summary 

In the present document we analyzed the behavior of the two-component plasma 
at T = 2 in the presence of one and two point adsorbing impurities, or in the 
presence of one and two adsorbing lines. 

For point impurities, as shown in Ref. [8], the partition function and corre- 
lations of the system can be expressed in terms of the same quantities for an 
unperturbed system, without impurities. Since at T = 2, exact results are avail- 
able for the partition function and correlations of the unperturbed system, we 
were able to obtain exact results for the partition function, the density profiles 
and correlation functions of the plasma with one or two point impurities. We 
also computed the electric potential created by one impurity, due to the charge 
it adsorbs and the polarization cloud that is formed around it. 

As a simple model for electrodes with adsorbing sites, we studied the prop- 
erties of the plasma with one or two parallel lines of absorbing impurities. The 
general formalism developed in Ref. [1] for the two-component plasma with an 
external potential is applied to this case. We obtained exact results for the par- 
tition function and density profiles of the plasma. In both cases we checked that 
various relations between the adsorbed charge and the excess grand potential 
are satisfied. 
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